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ABSTRACT 


We proposed and compared three methods (filter burnup, single-energy burnup, 
and burnup extremum analysis) to build a high-resolution neutronics model for ??*Pu 
production in high-flux reactors. The filter burnup and single-energy burnup methods 
have no theoretical approximation and can achieve a spectrum resolution of up to ~1 
eV, thereby constructing the importance curve and yield curve of the full energy range. 
The burnup extreme analysis method combines the importance and yield curves to 
consider the influence of irradiation time on production efficiency, thereby constructing 
extreme curves. The three curves, which quantify the transmutation rate of the nuclei 
in each energy region, are of physical significance because they have similar 
distributions. A high-resolution neutronics model for ??*Pu production was established 
based on these three curves, and its universality and feasibility were proven. The 
neutronics model can guide the neutron spectrum optimization and improve the yield 
of ?*Pu by up to 18.81%. The neutronics model revealed the law of nuclei 
transmutation in all energy regions with high spectrum resolution, thus providing 


theoretical support for high-flux reactor design and irradiation production of 7**Pu. 


Keywords: *°°Pu, Neutronics model, High-flux reactor, Spectrum resolution, Spectrum 


optimization. 


1. INTRODUCTION 


Plutonium-238 (?5Pu) is a radioactive isotope with a half-life of 87.7 years. Ри 
releases a-particles of 5.49 MeV, which are easily blocked with a range of 
approximately 20 um in an aluminum plate. Meanwhile, the a-decay of 2 Ри produces 
Uranium-234 (?*^U), which has a half-life of 2.455 х 10? years. Equation (1) gives the 
decay chain of ?5Pu. 

238 py (4) 87.7a 2341) (а) 2.455x105a 2307p, -— 206рр, (1) 
235Pu is an ideal heat source material without the need to consider the further decay of 
decay daughters [1]. 2 Ри has a density of 19.8 g/cm? (25 °C) and a heating power of 
0.57 W/g, and it has been widely used in radioisotope thermoelectric generators (RTGs) 
[2] and radioisotope heater units, such as the general-purpose heat source radioisotope 


thermoelectric generator (GPHS-RTG) in Galileo spacecraft [3] and the isotopic pulse 


cardiac pacemakers (IPCP) [4], as shown in Figure 1. 
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Figure 1. Structures of GPHS-RTG and IPCP 
High-flux reactors [5-6] have a stable and high neutron flux, facilitating ?*5Pu 
production. Two methods to produce ??*Pu are available: 1) in-reactor irradiation of 
Americium-241 (^! Am) [7] and 2) in-reactor irradiation of Neptunium-237 (77’Np) [8]. 
235Pu can be extracted from the decay products of Curium-242 (^? Cm) after the Ат 
target is irradiated in a reactor. Equation (2) provides the decay chain of this process. 


B^ 16h а 1634 
241 242 242 238 
gsAm(n, y )^9$5Am ———> 2 д С00---- Ри. (2) 


The purity of 23%и produced by this method is high, but this process releases strong 
gamma rays (59.3 eV), which deteriorate the radioactivity environment. There is almost 
no radioactivity problem during the production of ?*Pu by in-reactor irradiation of 
37Np. Therefore, the production of ?*Pu by in-reactor irradiation of 27Мр is a 
mainstream process [9]. Various nuclear reaction channels are coupled, and many new 
nuclides appear after the long-term irradiation of ?"Np targets, as shown in Figure 2. 


This study analyzed this production process. 
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Figure 2. Nuclide transmutation path of in-reactor irradiation of "Np 

The irradiation production of 2 Ри lacks a precise neutronics model, leading to а 
low transmutation rate of nuclei and high costs [10-11]. To improve production 
efficiency, the reactions occurring for different nuclides in the chain should be different 
during the irradiation period. For example, the target nuclide ?"Np should absorb as 
many neutrons as possible without any other reactions, the intermediate nuclide ?*Np 
should only undergo B-decay, and the produced ?*Pu should remain stable without any 
reaction occurring. It is difficult to regulate the nuclear reactions of the nuclide chains. 
The microscopic cross section of the nuclei 18 related to the energy spectrum [12]. 
Within some energy regions, the nuclide chain will have many ofthe required reactions, 
leading to a high efficiency of ?*Pu production, and vice versa. Therefore, neutron 
spectrum analysis and regulation can help increase the transmutation rates and reduce 
production costs. The optimization of irradiation in 2 Ри production can be 
theoretically divided into two topics: 1) determination of the optimal neutron spectrum, 


and 2) achievement of the optimal neutron spectrum. 
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Considerable research has been conducted to build neutronics models for the 
production of transuranic isotopes (such as ??*Pu and 22С#). Pan et al. [13] proposed a 
rapid diagnosis method for evaluating radiation schemes, which not only avoids tedious 
burnup calculations, but also helps provide direction for optimization. However, this 
method uses the initial nuclides in the target to represent all nuclides in the 
transmutation chain during the entire irradiation period, and only absorption and fission 
reactions are considered, while the other reactions (such as p-decay) cannot be 
considered. Pan et al. [14] defined key nuclides for neutron spectrum analysis, 
identified three energy regions that harm the transmutation of nuclides, and used 
filtering materials to reduce the neutron flux in these energy regions. However, this 
study did not consider the complete nuclei chain during irradiation. Hogle et al. [15-17] 
of the Oak Ridge National Laboratory (ORNL) presented a sensitivity curve detailing 
the production efficiency in each energy region. However, the analysis process was 
based on point-burnup calculations, resulting in a conclusion with no universality. 
Recently, Pan et al. [18] conducted a refined spectrum analysis of heavy nuclei 
synthesis in reactors, revealing the law of nuclear transmutation in all energy regions 
with a high spectrum resolution. However, this method has only been analyzed and 
verified for the synthesis of ?Cf and has not been applied to other transuranic isotopes. 
Therefore, refined neutronics models of ?*Pu production are lacking. 

We proposed and compared three methods for the high-resolution spectrum 
analysis of 23%и production in high-flux reactors. The transmutation rate of the nuclide 
chain in each energy region was quantified to build three relationship curves (referred 
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to as “importance curve,” “yield curve,” and “extreme curve") between the production 
efficiency and neutron spectrum. A high-resolution neutronics model for 238и 
production was established. The remainder of this paper is structured as follows: 


Section 2 introduces the spectrum analysis, Section 3 introduces the spectrum 


optimization, and Section 4 concludes the paper. 


2. SPECTRUM ANALYSIS 


2.1 Reactor and Target 

All the analyses were performed using a high-flux isotope reactor (HFIR) [19]. 
Rated at 100 MW and currently operating at 85 MW, the HFIR has a high steady-state 
neutron heat flux of 2.6 х 10? cm?:s'!. We used a pure ?"Np target, a 50 cm high and 
0.85 cm diameter rod placed in an irradiation channel at the neutron flux trap. The HFIR 
and target were modeled using the RMC code, a self-developed Monte Carlo code [20], 


as shown in Figure 3. 


reactor core neutron flux trap target 


Figure 3. Modeling diagram of the high-flux isotope reactor (HFIR) and the target 
2.2 Computational Method 
The yield of ?*Pu can be obtained by performing a Monte Carlo burnup 
calculation, which is a coupling of the Monte Carlo criticality and point-burnup 
calculations. The point-burnup equation describes the transmutation of nuclides over 
time when a target is irradiated in a high-flux reactor. For each nuclide in the burnup 


chain, the time-dependent point-burnup equation can be written as 


dn, eff 4 eff eff (1) 
- 28 n,—A, п, 
iz j 


where n; is the density of the i nuclide, A is the effective decay constant of the it 


nuclide, and Б, is the branching ratio for transmuting the i" nuclide to the j^ nuclide. 


A and b; can be calculated from the following formula: 


eff 
AU TA? 0. 
J 


b =b A eo, 4 AT 


EJ E 


(2) 


where 4; is the decay constant of the i^ nuclide, 0 is the neutron flux, and о; is the one- 
group cross-sections where the i" nuclide’s reaction generates the j nuclide. 
As shown in Equation (2), one-group cross-sections of the target are required; 


therefore, the Monte Carlo criticality calculation should be performed. 


1 
- E Es (3) 
(1-0-5)0 юН Еф 


where L is the leakage operator, C is the collision operator, S is the scattering operator, 
F is the fission operator, and keff is the effective multiplication factor. The physical 
parameters around the target, such as the neutron flux, fission reaction rate, and 
absorption reaction rate, can be obtained by solving Equation (3), which is used to 
determine the one-group cross-sections. 

The Monte Carlo burnup calculation requires one-group cross sections integrated 
according to the neutron spectrum over the entire energy range rather than the cross 
sections in each single energy region, which can only quantify the influence of the 
entire neutron spectrum on the production efficiency and not a closer analysis 
of each single energy region. Therefore, we cannot know which energy regions 
are favorable and which are unfavorable for ??*Pu production, and thus, cannot 
provide a more refined theory for optimization. 

We propose a filter burnup method and a single-energy burnup method for refined 
spectrum analysis, quantifying the transmutation rate of the nuclide chain in each 
energy region, and building relationship curves between the production 
efficiency and neutron spectrum. We also propose an extreme burnup analysis 
method to investigate the influence of the irradiation time on the relationship curves. 
Therefore, three relationship curves (importance, yield, and extreme) were 


obtained using the three analytical methods. 
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2.2.1 Filter Burnup Method 

The filter burnup method seems to perform subtraction on the neutron spectrum, 
i.e., reducing the neutron flux in a certain energy region to investigate the relationship 
between the change in grouped flux and the change in yield, quantifying the 
importance of the flux in a certain energy region on the transmutation rate of 
the nuclei. Because this is a sister paper to Reference [18], the detailed theory 


is not repeated and only the definition of importance is provided. 


AY,/Y (Y-Y)/Y (Y-Y)/Y 
Е E im (4) 
Аф, фф М -@ 


(12-21 
1 


where the subscript “i” is the energy region index, Y and Y; are the yields оЁ2 “Ри before 
and after flux reduction in the ith energy region, respectively, ф and 9 ' аге the grouped 
neutron flux before and after flux reduction, respectively, and M is the ratio of flux 
reduction. Previous analyses [18] have shown that the value of M (1/1, 1/2, 1/4, and 1/8) 
has little influence on the obtained importance curve; therefore, we take M — 1/8 directly. 

The entire energy range was divided into 238 regions [21]. The energy division for 
the filter burnup method has no limitations, thereby achieving a high spectrum 
resolution. The importance of each energy region throughout the 90-day irradiation 


period was calculated, building up the importance curve of 2 Ри production, as shown 


in Figure 4. The ten maxima and minima in the importance curve are listed in Table 1. 
1f 
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Figure 4. Importance curve of 2 “Ри production for the filter burnup method 
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Table 1. Maxima and minima in the importance curve 


Maximums Minimums 
Energy Regions (MeV) Values Energy Regions (MeV) Values 
[1.45х10%, 1.50x105] 1.91x10:? 11.85х107, 1.90х107| -2.34x10'!6 
14.50х107, 5.00х1077| 1.28х10:7 [2.77x10°, 2.87х10%] -1.33x10'!6 
15.00Х107, 5.50x107] 6.18x10:!6 [2.87x10 $, 2.97х10%] -7.466х1077 
[1.30x10°, 1.35x10°] 5.94x 10716 [5.00x107, 5.20x107] -4.64х 107" 
[1.40x10°, 1.4510] 5.18x107 [1.15x10^, 1.19104] -3.07х1017 
[1.50x10°, 1.59х10:| 3.18x10'!6 [1.00x10!, 1.28x10!] -2.59х1077 
[1.35x10°, 1.40x10°] 2.89х 1075 [8.19x109, 1.00x10!] -2.55x107!7 
14.00х107, 4.50x1077] 2.50х 1075 [6.43x109, 8.19x10?] -2.38х1077 
[1.20x10°, 1.50x10?] 2.21x107!6 11.70х107, 1.85х107| -2.29х1077 
[1.25x10°, 1.30x10°] 2.14х10:16 [1.85x109, 2.35x10?] -2.10х10:7 


The maximum values indicate that the neutrons in these energy regions are 
productive, that is, increasing the neutron flux in these energy regions will increase the 
yield of ?*Pu. Meanwhile, the minima indicate that the neutrons in these energy regions 
are harmful, that is, increasing the neutron flux in these energy regions will decrease 
the yield of **Pu. Therefore, we should increase the neutron flux in the energy regions 
with maximum values (referred to as the "positive energy region") and reduce the 
neutron flux in the energy regions with minimum values (referred to as the "negative 
energy region") to promote the production of ?*Pu. 

2.2.2 Single-energy Burnup Method 

The single-energy burnup method does addition on the neutron spectrum, 
assuming that the target was irradiated with a single-energy neutron source. We 
simulated the ?*Pu production efficiency with single-energy neutron sources, that is, 
the total neutron flux was the same, but all neutrons were in a certain energy region. 


The efficiency of a certain energy region was quantified by comparing the yields of 


2Pu with those of single-energy neutron sources. The full energy range was divided 
into 238 regions, and the yields of ?5Pu in each energy region throughout the 90-day 
irradiation were calculated, thereby constructing the yield curve for the single-energy 
burnup method, as shown in Figure 5. The ten maximum and minimum values in the 


yield curve are given in Table 2. 
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Figure 5. Yield curve of 2 Ри production for the single-energy burnup method 


Table 2. Maxima and minima in the yield curve 


Maximums Minimums 
Energy Regions (MeV) Values Energy Regions (MeV) Values 
[1.45x10 6, 1.50х10:1 1.17x10! [1.57x10!, 1.73x10!] 4.60x10^ 
(4.50х107, 5.00x107] 1.09x10! [1.73x10!, 2.00x10!] 4.60x104 
[1.30x10°, 1.35x10°] 1.01x10! [1.46x10!, 1.57x10!] 5.24x10^ 
[1.40x10°, 1.45x10°] 9.94x10° [1.38x10!, 1.46x10!] 6.04x10^ 
(5.00Х107, 5.50x1077] 9.84x10? [1.28x10!, 1.38x10!] 6.69x10^ 
[1.50х10°%, 1.59х10%] 9.10x10? [1.00x10!, 1.28x10!] 8.55x10^ 
[1.35x10 6, 1.40х10%] 8.89x10? [8.19x10°, 1.00x10!] 1.06x10? 
[1.25x10°, 1.30x10°] 8.29x10° [6.43x10°, 8.19x10?] 1.51x10? 
[3.73x10°, 4.00x 10°] 8.14x10? [4.80х10°, 6.43x10?] 2.97x10? 
(4.00х107, 4.50x107] 7.80x10? [4.30x10, 4.80x 10°] 3.99x10? 
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2.2.3 Burnup Extremum Analysis Method 

The importance curve obtained using the filter burnup method and the yield curve 
obtained using the single-energy burnup method exhibited similar variation trends in 
the resonance energy range, proving the physical nature of Figures 4 and 5. However, 
the specific values shown in Figures 4 and 5 are inconsistent. Eight of the ten maxima 
in Tables 1 and 2 correspond to the same energy regions, but only three of the ten 
minima, which is due to the different physical assumptions of the two analysis methods. 
The filter burnup method perturbs a particular energy spectrum to perform a 
perturbation analysis that describes the importance of a particular neutron spectrum, 
whereas the single-energy burnup method is not limited to a particular spectrum. Taking 
the high-energy regions as an example, the neutrons in these regions are almost useless 
for ?*Pu production. Therefore, the yield calculated by the single-energy burnup 
method is zero, and the importance calculated by the filter burnup method is also zero. 
However, the filter burnup method can also find the energy regions that are negative 
for °8Pu production under a particular spectral environment. Therefore, the importance 
can be negative numbers, which explains why only three of the ten minimum values in 
Tables 1 and 2 correspond to the same energy regions. 

The ?*7Np target [22] was irradiated in the neutron flux trap for 90 days using both 
the filter burnup and single-energy burnup methods. Significant differences in 238и 
production efficiency exist among different energy regions, resulting in different 
irradiation times required. For some energy regions, a 90-day irradiation exceeded the 
time required for the yield to peak, resulting in an increase first and then a decrease in 
the yield during the entire irradiation period (which is referred to as "excessive 
irradiation"), whereas in other energy regions, a 90-day irradiation cannot lead to peak 
production. We used the single-energy burnup method to calculate the yields of ?*Pu 
with 5-day and 10-day irradiation to investigate the influence of the irradiation time on 


the yield, as shown in Figure 6. 
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Figure 6. Comparison of curves with different irradiation times 
As shown in Figure 6, different irradiation times affected the curves. We proposed 
a burnup extreme analysis method to eliminate the influence of irradiation time on the 
importance curve (Figure 4) and yield curve (Figure 5), that is, the maximum derivative 
between the importance (or yield) and irradiation time was calculated and used to 
quantify the production efficiency in the 238 energy regions, building two extreme 


curves, as shown in Figure 7. 
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Figure 7. Two extreme curves 
It can be observed that the two curves have a similar distribution trend, regardless 
of the analysis method. Therefore, the importance curve in Figure 4, the yield curve in 
Figure 5, and the extreme curves in Figure 7 have physical significance and can jointly 


guide the spectrum optimization for 2и production. 
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2.3 Universality Testing 

All curves in Section 2.2 are calculated based on the HFIR, a thermal reactor. To 
test the universality of these curves, we performed the same calculations in another 
reactor with a fast spectrum [23-32], a high-flux lead-bismuth reactor (HFLBR). 
Detailed parameters of this reactor can be found in a previous study [13]. As previously 
discussed, the filter burnup method is based on a particular neutron spectrum, whereas 
the single-energy burnup method is not limited to a particular neutron spectrum. 
Therefore, the yield curve is independent of the reactor model, whereas the importance 
curve is related to it. The importance curve of the HFLBR for 23и production was 
calculated and compared to that of the HFIR, as shown in Figure 8. The ten maxima 


and minima in the importance curves are compared in Table 3. 
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Figure 8. Importance curves of the HFLBR and HFIR 
The two reactors exhibited different neutron spectra and flux levels. As shown in 
Figure 8, the two importance curves coincide well in the resonance-energy region and 
exhibit large deviations in the high- and low-energy regions, proving that the neutron 
spectrum influences the importance curve. However, as listed in Table 2, nine ofthe ten 
maxima and eight of the ten minima correspond to the same energy regions, proving 


that the importance curves are not completely dependent on the neutron spectrum and 
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are of physical significance, which can be used to determine the positive and negative 
energy regions. Therefore, the importance curve is universal and can be used for 
neutron spectrum optimization. 


Table 3. Maxima and minima in the two importance curves 


HFIR HFLBR 
Energy Regions Energy Regions 
(MeV) Values (MeV) Values 
[1.45x10°, 1.50x10°] 1.91х10:5 [1.45х10°%, 1.50x10°] 4.66107 
[4.50x107, 5.00107]  1.28x10'^ [4.50x10°7, 5.00107] 3.72х10:9 
[5.00x107, 5.50x107] 6.18x10°'° — [5.00x107, 5.50x107] (1.99х10:9 
[1.30x10°, 1.35x10 5] 5.94x10'!6  (1.40х109,1.45х109| 1.97x10:5 
[1.40x10°, 1.45x10 5] 5.18x10'6 — [1.30x10°, 1.35x10 5] 1.48х10:9 
Maximum 
[1.50x10°, 1.59x105] 3.18x10'6 — [1.50x10°, 1.59x10 5] 1.17х10:9 
[1.35x10°, 1.40x10°] 2.89x1075 [1.35х10°%, 1.40x10°]  1.01x10:5 
[4.00x107, 4.50x107] 2.50x1075 — [4.00x107, 4.50x107] — 7.45x10'!6 
[1.20x10?, 1.50x10?] 221x105 [1.25х10°%, 1.30x10°] 7.18х10:9 
[1.25x10°, 1.30x10 5$] 2.14x1076 . [3.73x10°, 4.00x10°°] 7.08х10:9 
[1.85x10°, 1.90x10?] -2.34x1076  [1.85x102, 1.90x10?] -7.63х10:6 
[2.77х10°%, 2.87х109| -1.33x10'6  [2.77x10/$, 2.87x10 5] -4.20х10:9 
[2.87x10°, 2.97x10$] -7.46x10"  [2.87x10$,2.97x10 5] -2.72х10:9 
[5.00x102, 5.20x10?] -4.64x10"  [1.28x10!, 1.38x10!] -2.34x10'!6 
11.15х107, 1.19x107] | -3.07x107 [1.70х10°, 1.85x10?] -1.18х10:6 
Minimum 
[1.00x10!, 1.28«10']  -2.59x107 [1.15х10*, 1.19х10*] -6.58х10"'” 
[8.19x109, 1.00x10!] -2.55х10!7  [2.67x10'6, 2.77x10'$] — -6.54x1077 
[6.43x109, 8.19x10?] | -2.38x10"  [1.00x10!, 1.28x10!] -5.02х1077 
[1.70x10°, 1.85х107| -2.29x10"" — [8.19x10*, 1.00x10!] -476х1077 
[1.85x109, 2.35x10°] -2.10х1077 — [6.43x109, 8.19x10°] -447x10" 
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3. VERIFICATION AND APPLICATION 


3.1 Neutronics Model Verification 

We obtained the importance curve (Figure 4, marked as '/;") by the filter burnup 
method, the yield curve (Figure 5, marked as *Y;") by the single-energy burnup method, 
and the two extremes curves by the burnup extremum analysis method (Figure 7, 
marked as “Е; (/;)” and “Е; (Yi)"). These curves were used to build a high-resolution 
neutronics model for 238и production in high-flux reactors, which could be used to 
guide neutron spectrum optimization to improve production efficiency. 

To demonstrate that these curves can be used to guide the neutron spectrum 
optimization for ?*Pu production in high-flux reactors, we constructed a variety of 
irradiation schemes by dispersing nuclides into the target to obtain many different 
neutron spectra. We determined whether the yields of these irradiation schemes were 
positively correlated with the corresponding total spectral efficiency. The total spectrum 


efficiency (marked as “7”) of an irradiation scheme was calculated as follows: 
T-5 4X, (5) 


where N is the total number of divided energy regions, 0, is the neutron flux in the i” 
energy region, X; represents the efficiency of the i™ energy region, preferably Ji, Y;, E; 
(1), or Е; (Yi) where T corresponds to Ті, 7, Тз, and Та. 

The nuclides we selected included !°’Ag, !?Eu, !°’Gd, ?Tm, !é! Dy, 3 Ey, 5$, 
108рд, 140Ва, 40Ar “Li, “Ni, Ti, !5!Eu, “Не, "Та, Sm, and '8°W. The number of 
nuclides added could be 107, 107, 107, and 10? (іп (10*4#/cm’)). A total of 56 
irradiation schemes were constructed. The correlation coefficients between AY and AT 


were calculated as follows: 


Cov(AY, AT,) 


AY, AT,) = ===, 
" : JVar[AY WVar[AT, | 


п-1,2,3,4, (6) 
where AY is the variation of 23%и yield, AT is the variation of the total spectrum 
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efficiency, Cov(x,y) is the covariance of x and y, and Var[x] is the variance of x. The 
variation in the total spectrum efficiency T of these radiation schemes and the 


corresponding variation in 238и yield are listed in Table 4. 


Table 4. Variation in T and corresponding variation in “Ри yield 


Nuclides Amount AY(%) ATi АТ» АТ, АТ: 

Ag 104 0.05 8.45x10? -3.83x1013 4.78x10? 7.92x10? 
10° 0.04 5.82x10? -1.82х10'3 3.60x10? 3.38x10? 
йл 10^ -1.50 -4.96х102 -1,33х109 -4.24х102 -1.95х10" 
10° 0.21 1.43x10? -2.31х10'4 -24/7x10? -2.57х10!0 
10+ -11.62 -3.31х10"! -7.70x10!5 -2.79x10! — -1.23x10? 

1500 
10° -2.48 -8.80х102 -1.93х10'% -7.61х102 -3.25х10" 
T 10+ -0.20 5.40x10+ -9.66х10'% 1.91x10^ -6.89x10? 
10° -0.07 6.32x10? -5.90x105? 2.68x10? 7.37х108 
ару 10^ -0.11 -1.81x10? -1.47x10'4 -2.23x104 . -1.20x10!° 
10° -0.16 3.14x10? -1.37х10'4 -5.80x10^ -1.04х10!0 
brin 10^ -0.07 5.08x10? -1.14x104 2.22x10? -2.26x10? 
10° -0.07 3.67x 10° -2.80х10'3 1.60x10? 2.15x105 
Be 10+ -0.22 -5.32х102 -1.50x10!5 -5.04х102 — -229x10!! 
10° -0.16 -2.14x10? -2.46x10'4 -4.57x107  -3.07x10!? 
тэр 10+ -0.08 4.08x 10° -3.63x10% 1.16x10? -1.69x10? 
10° 0.05 4.02x104 -4.52x 10! 5.64«104 -1.09x10? 
10? -0.39 -4.53х102 -1.90х10' -1.84х102 -1.71х10" 
lias 10? 0.32 -3.05x10? -4.52x 1014 4.27x104 -3.10x10!? 
10^ -0.06 4.94x10? -1.13x104 2.72x10° -3.80x10? 
10° -0.05 1.54x10? -4.37x1013 1.57x10? 6.16x105 
10? -0.68 9.07x10? 1.121033 4.87x10? 9.10x10? 
зу 10? -0.11 9.96x 104 -7.68x1013 1.13x10° -1.68x10? 
10+ -0.02 1.82x10? 1.30x10! 2.73х10% 3.75x10? 
10° 0.01 1.68107 -1.04х10' 1.95x10? -4.00x10? 
10? -0.70 9.71x104 2.69x 10% 1.00x10? 8.71x10? 
m 103 -0.08 -4.37x10? -1.04х10' 4.62x104 -6.63x10? 
10+ -0.10 8.55x 103 -4.90x 10? 4.14x10? 2.04x10? 
10° -0.11 1.51x104 -4.18х10" -1.09x104 -3.68x10? 
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10? -0.70 1.25x10? 2.92x10! 3.56x10? 8.88x10? 


Bs 10° -0.14 6.81x10? -1.08x10'4 1.20x104 -1.02x101? 
10^ -0.08 -2.32x104 -2.80*1013 6.53x104 2.87x10? 
10° -0.11 4.52х10* -1.07x10" 1.81x107 -7.04x10? 
10? -0.61 3.98x10? 8.77х10'3 4.48x10? 1.19x101? 

i 10? 0.01 2.94x10? -2.84x10!! 3.28x10? 6.06x10? 
104 -0.06 7.71x10° -9.80х 1012 8.98x10^ 3.24x108 
10° -0.05 6.60x10? -3.74x105 1.20x10? -1.09x10? 
10? -1.34 -4.43х102 -1.15x10!° -3.48х102 -1.62х10" 

7 10? -0.24 1.80x10? -1.89x10" -7.86x10^ -1.66х10'% 

-= 104 -5.46 -1,75х101 -4.13х109 -1.34x107 -6.03x10!! 
10° -0.69 -2.38x10? -6.42х10'* -1.55x10? -8.29х1010 

102 -0.66%  1.07x10? 1.5910" 1.76x103 4.09x 10? 

siis 10° -0.06 -2.15x10? -1.25x10'4 8.74x 10" -6.80x10? 
104 -0.13 -1.99x 104 -1.12x10" 6.50x10^ -8.23x10? 

10° -0.18 -1.90x103 -9.50x 1013 -9.98x10^ . -1.12x10!? 

10? -1.33 -2.73х102 -1.21x10!5 -1.80х102 -1.34х10" 

"A 10? -0.07 6.75х10* -1.96х10'* 1.46x107 -1.99x10!° 
10^ -0.02 5.25x10? -6.28х1013 2.54x10? -2.08x10? 
10° -0.14 -2.00x10? -5.80*1013 2.16x10? 2.08x10? 
10? -0.39 -4.53x10? -1.90х10' -1.84х102 -1.71х10" 

T 10° 0.32 -3.05x10? -4.52x10!4 4.27x10* -3.10x10!° 
104 -0.03 9.17x10? -1.39x10!4 1.12x10? -1.05x10!? 
10° -0.10 3.78х10* -5.79х10' 1.74x103 -1.98x10? 
10? -0.88 -3.10х102 -9.95x10'4 -1.54х102  -1.10x10!! 

if 10? 0.45 1.12х102 -1.65x10" 4.10x103 -1.02x101? 
104 -0.03 1.03x10? -2.71х108 3.16x10? -1.3х10° 
10° -0.08 4.75x10? 2.90« 10% 2.10x10? 2.67x10? 


correlation coefficients 


between AY and AT 


0.949 -0.452 0.933 0.139 


As can be observed from Table 4, as the filter burnup method is based on a 
particular neutron spectrum, AY is positively correlated with AT; and A75, with 


correlation coefficients larger than 0.9. The single burnup method is not based on a 
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particular neutron spectrum, and hence AY is not positively correlated with A7? and А74. 
Therefore, the importance curve based on the filter burnup method and the extreme 
curve based on the importance curve can be used to guide the optimization of a 
particular irradiation scheme, that is, to further increase the yield оГ2 Ри of a particular 
irradiation scheme. The yield curve based on the single-energy burnup method and the 
extreme curve based on the yield curve cannot be used to guide the optimization of a 
particular irradiation scheme. However, the single-energy burnup method can better 
determine the positive- and negative-energy regions, particularly when the neutron flux 
in the energy region accounts for a small amount of the total neutron flux. Therefore, 
the single-energy burnup method can amplify the difference in importance of the filter 
burnup method, which is not sufficiently explicit in some energy regions. Consequently, 
the two methods complement each other. 

3.2 Neutronics Model Application 

As shown by the importance, yield, and extreme curves, neutrons in low-energy 
regions are conducive to the production of 2 Ри, Therefore, we slowed the neutrons in 
the target by dispersing 'H in the target and constructing thermalized neutron energy 
spectra. The amount of hydrogen added started from 0.0 (1074#/cm?), increased by 0.2 
(1074#/cm?) each time, and ended at 3.0 (1074#/cm?), resulting in a total of 16 irradiation 
schemes (including the original scheme). The results of AY(%), Tı and 73 are presented 
in Figure 9, drawn as a dot plot. 

As shown in Figure 9, the yield of ?*Pu can be increased by regulating the neutron 
spectrum around the target. In the above 16 irradiation schemes, the rate of increase in 
the yield of ?5Pu was the highest when the amount of ІН added was 2.6 (1029 си), 
reaching 18.81%. Meanwhile, the rate of increase was positively correlated with the 
total spectrum efficiency, implying that the yield of ??*Pu can be further improved by a 
more refined neutron spectrum optimization. Therefore, a high-resolution neutronics 


model based on the importance curve, yield curve, and two extreme curves can guide 


18 
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the neutron spectrum optimization to improve the production efficiency of 238Ри. 
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Figure 9. Results of AY(%), T1, and 73 
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4. CONCLUSION 


?38Py is an ideal heat source material widely used in radioisotope heater units. The 
production of **8Pu by in-reactor irradiation of ?"Np is a mainstream process that 
involves complex nuclear transmutation. The irradiation production of ?*Pu lacks a 
precise neutronics model, which leads to a low transmutation rate of the nuclei and high 
costs. We proposed and compared three methods (filter burnup, single-energy burnup, 
and burnup extremum analysis) to build a high-resolution neutronics model for ??*Pu 
production in high-flux reactors. 

The filter burnup and single-energy burnup methods have no theoretical 
approximation and can achieve a spectrum resolution of up to ~1 eV, thereby 
constructing the importance and yield curves of the full energy range. The filter burnup 
method is based on a particular neutron spectrum, whereas the single-energy burnup 
method is not limited to a particular neutron spectrum. The burnup extreme analysis 
method combines the importance and yield curves to consider the influence of 
irradiation time on production efficiency, thereby constructing extreme curves. The 
three curves, which quantify the transmutation rate of the nuclei in each energy region, 
are of physical significance because they have similar distributions. A high-resolution 
neutronics model for ??*Pu production was established based on these three curves, and 
its universality and feasibility were proven. 

The neutronics model can guide the neutron spectrum optimization and improve 
the yield of ?*Pu by up to 18.81%. The neutronics model revealed the law of nuclei 
transmutation in all energy regions with high spectrum resolution, providing theoretical 
support for high-flux reactor design and irradiation production of 2 Ри, Moreover, the 
flowchart of the three methods for spectrum analysis is universal and can be used for 
the spectrum analysis of other scarce isotopes. In the future, we will conduct research 
on the following three aspects: (1) target design from an engineering perspective, (2) 
spectrum optimization based on the proposed neutronics model, and (3) irradiation 


channel construction for ?*5Pu production in a high-flux reactor. 
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